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1. Background 

1.1 Cosmology with weak leasing surveys 

Weak gravitational leasing is the coherent distortion in the shapes of distant galaxies due to 
the deflection of light rays by intervening mass distributions. Measurements of the effect on large 
scales is termed “cosmic shear” and has emerged as a powerful probe of late-time cosmology over 
the last 15 years (see e.g. Heymans et al. 2013 for recent results from the CFHTLenS survey). 
Since gravitational leasing is sensitive to the total (i.e. dark plus baryonic) matter content of the 
Universe, it has great potential as a very robust cosmological probe, to a large degree insensitive to 
the complications of galaxy formation and galaxy bias. One of the most promising aspects of weak 
leasing measurements is their combination with redshift information: such measurements are then 
a sensitive probe of both the geometry of the Universe and of the evolution of structure over the 
course of cosmic time. In turn, these latter effects are dependent on the nature of the dominant dark 
energy component in the Universe and/or on modifications to the theory of General Relativity on 
large scales. 

The observed distortions in the shapes of distant galaxies yields an estimate of the leasing 
shear field, 7. Since gravity is a potential theory, the shear at angular position d can be related to a 
leasing potential (i/r) as 

Yij{e) = (^didj- ¥{d), ( 1 - 1 ) 

where 5,- = r{5jj — rirjVi) is a dimensionless, transverse differential operator, and 5^ = 5i5^ is 
the transverse Laplacian. The indices {i,j) each take the values (1,2). In equation (1.1) we have 
assumed a flat sky which is an excellent approximation for the scales of interest (i.e. from ~ 100 
/j^^kpc to ~ 100 /i^^Mpc). The leasing potential can in turn be related to the 3-d gravitational 
potential, 4>(r) by (e.g. Kaiser 1998; Hu 2000) 



where r is the comoving distance to the sources. 

In the limit of weak leasing (7 << 1), and in the absence of intrinsic correlations in galaxy 
ellipticities, one can form an unbiased estimator for the shear field at a given sky position from 
the average of the observed galaxy ellipticities at that position. Since the shear field induced by 
large scale strucutre is small (typically ~ a few %) compared to the intrinsic dispersion in galaxy 
ellipticities (Oe ~ 0.3) one needs to average over many background galaxies in order to obtain a 
precise measurement. Measurements of cosmic shear can be affected by a number of instrumental 
and astrophysical systematic effects. The primary insturmental effect of concern is anisotropies in 
the point spread function (or beam) of the telescope which can mimic a cosmic shear signal. On 
the astrophysical side, alignments in the intrinsic shapes of galaxies can also mimic a cosmic shear 
signal. Consequently a great deal of effort is currently focused on careful instrument design for the 
next generation of weak leasing surveys as well as theoretical modelling and numerical simulations 
of intrinsic alignment effects. 

Observationally, to date the field of weak leasing has largely been the preserve of optical 
surveys due to the much larger number densities of background galaxies achieved in such surveys. 
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However, this will ehange with the advent of the Square Kilometre Array whieh will reach number 
densities of well-detected and well-resolved galaxies of up to ~ 5 galaxies arcmin^^ over several 
thousand deg^ in Phase 1, and ~ 10 galaxies arcmin^^ over 371 steradians in Phase 2. In addition, 
as described in this chapter, the radio offers truly unique approaches to measuring weak lensing 
that are (i) not available to optical surveys and (ii) potentially extremely powerful in minimizing 
the most worrying systematic effects in weak lensing cosmology. 

SKA surveys will also extend the reach of weak lensing beyond that of the optical mega¬ 
surveys of LSST and Euclid. First the typical redshifts probed by SKA surveys will go beyond that 
of optical surveys. Secondly the SKA offers the prospects of measuring the weak lensing distortion 
in 21cm HI intensity maps. Such an intensity mapping lensing survey will not only extend to the 
high-redshifts inaccessible to optical surveys but will also include precise redshift information to 
accompany the distortion signal. This signal will fill the gap between traditional galaxy lensing 
signal (where the sources are typically located at z ~ 1) and the CMB lensing signal originating 
at z = 3000. It thus holds the promise of yielding very precise cosmological constraints on the 
evolution of structure during early times where structures are better described with linear physics 
and can provide unique insight into so-called “early dark energy” models which exhibit observable 
signatures at early times. 

1.2 Radio weak lensing studies to date 

To date, the application of weak lensing analyses to radio data has been rather limited due to 
the relatively low number density of background galaxies achieved in large scale radio surveys. 
The only major weak lensing analysis of a radio survey was performed in Chang et al. (2004) who 
detected a cosmic shear signal in the Faint Images of the Radio Sky at Twenty cms (FIRST) survey 
(Becker et al. 1995), conducted with the Very Farge Array (VFA). Although the number density 
of sources in FIRST was low by optical standards (FIRST contains ~ 90 sources deg^^ compared 
with typically ~ 10 sources arcmin^^ in deep optical lensing surveys), a detection of cosmic shear 
on large scales was achieved by virtue of the large survey area covered (~ 10,000 deg^). 

The analysis of Chang et al. (2004) made use of the shapelets method (Chang & Refregier 
2002; Refregier & Bacon 2003) to measure the shapes of the FIRST sources directly from the uv 
interferometric data and included a thorough treatment of the possible systematics that can affect 
radio-based lensing shear estimates in particular. 

More recently, Patel et al. (2010) revisited the idea of radio-based lensing shear estimation 
through a re-analysis of the combined VFA -i- MERFIN observations of the Hubble Deep Field 
North (Muxlow et al. 2005). This work presented an implementation of image-based shapelets for 
shear estimation and also highlighted the advantages of cross-correlating radio-based and optical- 
based shear estimates obtained over the same area of sky. 

1.3 The path to SKA weak lensing 

A great deal of algorithm development and new analysis techniques are required in order to 
develop the field of radio weak lensing beyond the initial studies mentioned above. Key areas 
where work is required are the development of shape and/or shear estimation techniques that are 
well suited to radio interferometer datasets, the development of cross-correlation techniques for 
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combining optical and radio-based shear estimates and the development and demonstration of novel 
radio-specific lensing approaches making use of polarization and rotational velocity information. 
A number of SKA pathfinder and precursor telescopes lend themselves naturally to developing 
several of these areas. 

A key pathfinder telescope for demonstrating the long-baseline high resolution observations 
suited for weak lensing is the e-MERLIN interferometer based in the UK. A number of e-MERLIN 
legacy programs are well-suited to pathfinding the techniques for SKA weak lensing. Of particu¬ 
lar note are the e-MERGE* and SuperCEASS^ projects. e-MERGE is a multi-tiered project that 
(amongst other things) will image the GOODS-N field to 0.5 /iJy rms in the central 100 arcmin^ and 
to 1 /xJy rms in the surrounding 800 arcmin^. The SuperCEASS project complements e-MERGE. 
Its primary science driver is to detect the weak lensing effects of a supercluster of galaxies lo¬ 
cated at z = 0.2. SuperCEASS will image a 1.75 deg^ field to 4 /rJy rms. Both e-MERGE and 
SuperCEASS will perform their primary observations at 1.4 GHz with a resolution of 200 mas. 
Both will therefore act as training experiments for demonstrating shape measurement and shear 
extraction algorithms on high-resolution radio data, on the path to the SKA. 

A further opportunity exists to refine and test radio lensing analyses with the upgraded JVEA. 
A series of sky surveys (VEASS^) are currently being discussed for the JVEA and these will po¬ 
tentially include a deep fields component. If conducted with an array configuration that includes 
long baselines, these could be extremely powerful for pathfinding weak lensing techniques in the 
radio over larger survey areas, approaching 10 deg^ (Brown et al. 2013). 

Other opportunities for demonstrating weak lensing techniques in the run up to the SKA 
include EOEAR surveys^ (when the international baselines are included) and the CHIEES^ and 
CHIEES con Pol'^ surveys on the JVEA. The EOEAR surveys will be critical for testing weak 
lensing techniques at low frequencies while the CHIEES and CHIEES con Pol surveys will prove 
useful for testing novel ideas for radio lensing including the use of polarization information and HI 
rotational velocity measurements (see Section 2.4). 

2. Cosmic shear with radio continuum surveys 

2.1 Accessing the largest scales and the highest redshifts 

One of the potentially most exciting aspects of weak lensing surveys with the SKA will be the 
extra reach in terms of accessing the largest scales in the Universe and going beyond the reach of 
other weak lensing surveys in terms of the redshifts that are probed. We illustrate this in Pigs. 1-3 
where the redshift distributions of weak lensing source galaxies and the corresponding forecasted 
errors on a set of tomographic shear power spectra are presented. We present these forecasts for a 
2-year continuum survey with the SKA at three different stages of development - an early phase 
of SKAl comprising 50% of the envisaged sensitivity targeting a survey area of 1000 deg^; the 

'http://www.e-merlin.ac.uk/legacy/projects/emerge.html 

^http://www.e-merlin.ac.uk/legacy/projects/superclass.html 

^https://science.nrao.edu/science/surveys/vlass 

^http://lofar.strw.leidenuniv.nl 

^http://www.mpia-hd.mpg.de/homes/kreckel/CHILES/index.html 
®http://www.aoc.nrao.edu/~chales/chilesconpol/ 
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Table 1: Observational parameters used to produce the power spectrum and cosmological parameter fore¬ 
casts in Figs. 1^. 


Survey 

A,ky (deg2) 

tigai (arcmin ^) 

Zm 

SKAl-early 

1000 

3.0 

1.0 

VST-KiDS 

1500 

7.5 

0.6 

SKAl 

5000 

2.7 

1.0 

DES 

5000 

6.0 

0.6 

SKA2 

30940 

10 

1.6 

Euclid 

15000 

30 

0.9 


full SKAl surveying 5000 deg^, and a SKA2 survey covering 371 steradians. To generate these 
forecasts, we have adopted the performance specifications for SKAl outlined in Braun (2013). 
In particular, we have modeled the performance of Band 2 of the SKA-Mid facility as under the 
current baseline design, this telescope and frequency band combination provides the most powerful 
survey speed for the high angular resolution observations required for weak lensing. We have 
further assumed an object detection threshold of S/N > 10 and an angular resolution requirement 
of ©res = 0.5 arcsec. To model the redshift and flux dependence of the source population, we have 
made use of the SKA Design Studies (SKADS) simulations of Wilman et al. (2008), updated to 
match the galaxy number counts observed in the deepest radio surveys performed to date (Muxlow 
et al. 2005; Morrison et al. 2010; Schinnerer et al 2010). We have also assumed an RMS dispersion 
in intrinsic galaxy ellipticities of /nns = 0.3. 

For comparison on these plots, we also show the corresponding forecasts for optical weak 
lensing surveys that will be conducted over similar survey areas and on comparable timescales. 
Specifically, we consider the VST-KiDS (de Jong et al. 2013), the Dark Energy Survey^ and the 
Euclid satellite mission (Eaureijs et al. 2011). The observational parameters adopted to produce 
these forecasts are summarized in Table 1. In all cases, the radio surveys extend to higher redshift 
than the corresponding optical probes. They thus hold the potential to probe the power spectrum at 
higher redshift providing a more sensitive lever arm with which to constrain the growth of structure 
over cosmic time. 

Eig. 4 presents forecasted constraints on the matter density (D^) and matter power spectrum 
normalization (ag) cosmological parameters for the same envisaged SKA surveys as were adopted 
to generate Eigs. 1-3. Note that these forecasts are presented for the case of the standard 6- 
parameter ACDM model and no prior information is assumed - that is the projected constraints 
are coming solely from the envisaged SKA weak lensing survey. 

To generate the constraints, we have computed a simple shear power spectrum covariance 
matrix from Takada & Jain (2004), and we use the COSMOSIS cosmological parameter estima¬ 
tion code (Zuntz et al. 2014) to compute power spectra, parameter constraints, and marginalised 
contours. Note that no systematic errors are included in the analysis; errors are purely statistical. 
However, we have attempted to take into account anticipated knowledge and uncertainties regard- 

^http://www.darkenergysurvey.org 
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Figure 1: Left panel: The redshift distribution of source galaxies for a 1000 deg^ weak lensing survey 
requiring 2 years observing time on the SKA 1-early facility. Also shown is the redshift distribution for the 
1500 deg^ VST-KiDS optical lensing survey. The n{z) extends to higher redshifts in the radio survey and 
probes a greater range of cosmic history. Right panel: The corresponding constraints on a 5-bin tomographic 
power spectrum analysis. For both experiments, we assumed an RMS dispersion in ellipticity measurements 
of 7rms = 0.3 and the tomographic bins have been chosen such that the bins are populated with equal numbers 
of galaxies. Note how the radio survey extends to higher redshifts where the lensing signal is stronger and 
therefore easier to measure. Open triangles denote Iff upper limits on a bandpower. Note that only the auto 
power spectra in each bin are displayed though much cosmological information will also be encoded in the 
cross-correlation spectra between the different z-bins. 



Figure 2: As Fig. 1 but for a 5000 deg^ weak lensing survey requiring 2 years observing time on the 
full SKAl facility. Also shown for comparison are the n{z) distribution and forecasted power spectrum 
constraints for the 5000 deg^ Dark Energy Survey. 


ing photometric and spectroscopic redshift estimates for the background galaxy population. For 
SKAl-early, we have assumed that we have no spectroscopic redshift information and that we have 
photo-z estimates from overlapping optical surveys with errors = 0.05(1 +z) up to a limiting 
redshift of 1.5. To model the much larger uncertainties expected for the high-z radio galaxies, we 
adopt ff^ = 0.3(1 +z) so that a z = 2 galaxy has a redshift uncertainty of ± ~ 1. For SKAl, we 
additionally assume that we will have spectroscopic redshifts from overlapping HI observations 
for 15% of the z < 0.6 population. Finally for SKA2, we assume we have spectroscopic redshifts 
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Figure 3: As Fig. 1 but for a 37Z steradian weak lensing survey requiring 2 years observing time on SKA2. 
Also shown for comparison are the n{z) distribution and forecasted power spectmm constraints for the 15000 
deg^ Euclid satellite mission. 

for 50% of the z < 2.0 population. The forecasts presented in Fig. 4 account for these redshift 
uncertainties. 

We see from Fig. 4 that even the SKAl-early survey targeting the smallest sky area can provide 
competitive constraints on cosmological parameters — the forecasted constraints for the 1000 deg^ 
SKAl-early survey are a factor of ~5 better than the tomographic weak lensing analysis of the 
current state-of-the-art CFHTLenS data (Heymans et al. 2013). We also see large improvements in 
the constraints obtainable with each subsequent stage of the SKA — the constraints obtainable with 
SKAl are broadly comparable with the KiDS and DBS optical surveys while SKA2 is competitive 
with Euclid. 

Fig. 4 also demonstrates that our nominal choice of target survey areas for the three stages of 
the SKA are, broadly speaking, optimal choices from the point of view of constraining these cos¬ 
mological parameters — for SKAl-early the 1000 deg^ survey provides the strongest constraints, 
for SKAl the 5000 deg^ survey performs the best while for SKA2, the 371 steradian survey provides 
the best constraints. 

2.2 The promise of radio observations to suppress weak lensing systematics 

Optical and radio surveys, such as Euclid and/or ESST and the SKA, have a particularly useful 
synergy in reducing and quantifying the impact of systematic effects which may dominate each 
survey alone on some scales. By cross-correlating the shear estimators from one of these surveys 
with those of the other, several systematic errors are mitigated. 

We can see this by writing the contributions to an optical (o) or radio (r) shear estimator: 



( 2 . 1 ) 

( 2 . 2 ) 


where /grav is the gravitational shear we are seeking, is the intrinsic ellipticity of the object. 


and 7sys are systematic errors induced by the telescope. If we correlate optical shears with optical 
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Figure 4; Forecasted constraints in the ds — Q.m parameter space for 2-year continuum surveys using SKA- 
Mid/Band 2 with SKAl-early {upper left panel), SKAl {upper right panel) and SKA2 {lower panel) perfor¬ 
mance parameters. For each of these cases, we present the constraints for survey areas of 1000 deg^ (blue) 
5000 deg^(red) and 2>n steradians (green). 


shears, or radio shears with radio shears, we obtain terms like 


(77) ~ ( 7 grav 7 grav) + ( 7 grav 7 int) + ( 7 int 7 int) + ( 7 sys 7 sys)) 


(2.3) 


where the first term is the gravitational signal we seek, the second term is the GI intrinsic alignment 
(Hirata & Seljak 2004), the third term is the II intrinsic alignment (e.g. Heavens et al. 2000), and the 
final term is the contribution from systematics. All of these terms could be similar size on certain 
scales, which is damaging to cosmological constraints. On the other hand, if we cross-correlate the 
optical shears with radio shears, we obtain 


(yWyf'-)) = (7g„v7grav) + (7grav7S) + (7grav72) + (7^7^) + (7^7^) 


(2.4) 


The second and third terms are the GI alignment (Hirata & Seljak 2004), which still survives. 
However, the fourth term involves the correlation between optical and radio shapes, which will be 
less than that between one frequency alone as the emission mechanisms are different (c.f. Patel 
et al. 2010 where no correlation at zero lag was found). This term is therefore reduced. Most 
importantly, the fifth term involving systematics is expected to be zero, as the systematics in these 
two telescopes, which are of completely different design and function, are not expected to be 
correlated at all. We are therefore able to remove the dangerous systematics correlation from our 
shear analysis - and to gain an estimate of its magnitude in the autocorrelation case. 
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Table 2: Requirements on multiplicative and additive biases on ellipticity measurement for proposed SKA 
weak lensing surveys to be dominated by statistical rather than systematics uncertainties. Q is a global 
“quality factor” which we calculate from m and c following Voigt et al. (2010). 


Experiment 

^sky 

^gal 

Zm 

m < 

c < 

Q> 

SKA 1-early 

1000 

3.0 

1.0 

0.014 

0.0012 

62 

SKA 1-early 

5000 

1.2 

0.8 

0.012 

0.0011 

79 

SKA 1-early 

30940 

0.35 

0.5 

0.011 

0.0011 

80 

SKAl 

1000 

6.1 

1.2 

0.0090 

0.00095 

103 

SKAl 

5000 

2.7 

1.0 

0.0067 

0.00082 

140 

SKAl 

30940 

0.9 

0.7 

0.0058 

0.00076 

164 

SKA2 

1000 

37 

1.6 

0.0031 

0.00055 

318 

SKA2 

5000 

23 

1.4 

0.0019 

0.00043 

523 

SKA2 

30940 

10 

1.3 

0.0012 

0.00035 

825 


A further systematic that is potentially problematic for weak lensing surveys are color-gradients 
(Voigt et al. 2012). However, galaxies typically have smooth spectral dependence at radio fre¬ 
quencies. This, combined with the ability to measure the spectral dependence of the beam very 
accurately, means that radio weak lensing surveys will be insensitive to this systematic. 

2.3 Shear measurement 

In order to extract the weak lensing science from the SKA, we will need to make highly accu¬ 
rate measurements of the shearing of background sources in real, noisy data. In optical experiments, 
this shear measurement process has typically consisted of measuring the ellipticity of individual 
galaxies and combining these measurements to form an estimate of the cosmic shear. A wide vari¬ 
ety of galaxy shape measurement algorithms have been developed for this task and tested through 
the STEP and GREAT programmes (see e.g. Mandelbaum et al. 2014 and references therein) and 
experience with real data, honing the techniques and focusing research on areas which require im¬ 
provement. As a consequence, shape measurement from optical data is a well-established field and 
techniques already developed are sufficiently advanced that shape measurement-induced system¬ 
atics are likely to be sub-dominant to statistical errors in current and near future optical lensing 
surveys. These shape measurement systematics are often parameterised in terms of the additive 
bias c and multiplicative bias m on the ellipticity recovered from an input source with known 
true ellipticity 

= (l+my™^ + c. (2.5) 

A formalism is provided in Amara & Refregier (2008) for calculating the m and c necessary for 
systematic errors to be sub-dominant, with requirements for a selection of SKA surveys shown in 
Table 2. However, these algorithms have almost solely been motivated by optical and NIR data, 
meaning their robustness to issues peculiar to radio data, in particular potential systematics from 
the non-linear deconvolution process necessary for imaging (and removal of sidelobe aitifacts), is 
unclear. Of the two radio weak lensing studies so far performed, both measure shear by modeling 
data using the orthonormal shapelets basis set (Refregier & Bacon 2003). Patel et al. (2010) use 
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shapelets to model images reeonstrueted using the CLEAN algorithm, while Chang et al. (2004) 
take advantage of the faet that shapelets remain loealised analytie funetions under Fourier trans¬ 
formation to direetly model the visibility plane and were able to make a 3.6a deteetion of eosmie 
shear. More reeently, Patel et al. (2013) have applied image plane shapelets to simulations, allow¬ 
ing them to quantify the effieaey of the method, aehieving values m = 0.176 and c = 0.006 for a 
simulated e-MERLIN observation. For SKAl, Patel et al. in this volume aehieve m = 0.28 and 
c = 0.001 for image plane shapelets, eomparable to requirements of m < 0.0054, c < 0.0073 for 
the SKAl 5000 deg^ survey, as well as providing an initial eomparison of the relative performanee 
of visibility plane shapelets. 

A number of open problems still remain. Although initial steps have been taken, only a small 
part of the spaee of possible radio shear measurement teehniques has so far been explored. Future 
work may be expeeted to take plaee with two main themes: 

• Investigation of the ability of image reeonstruetion algorithms to produee images with the 
level of fidelity neeessary for weak lensing eosmology and subsequent adaptation of known 
methods from optieal/NIR weak lensing to radio data. High-fidelity image reeonstruetion is 
an aetive researeh topie in its own right, with many extensions and alternatives to traditional 
CFEAN and Maximum Entropy methods under development (e.g. Sutter et al. 2014; Carrillo 
et al. 2014 and referenees therein). 

• Investigation of teehniques whieh measure shear direetly from interferometer visibilities. 
Sueh teehniques have the advantage of avoiding systematies introdueed by the non-linear 
deeonvolution neeessary for the imaging proeess and a visibility plane teehnique (shapelets 
in Chang et al. 2004) has been the only teehnique so far to sueeessfully deteet weak lensing 
in radio data. However, for the datasets produeed by the SKA, eomputational ehallenges 
for visibility plane teehniques will be great. Simultaneously fitting ~5 parameter models to 
large numbers of sourees over large numbers of visibilities is likely to be unfeasible, though 
averaging and novel methods sueh as uv-staeking (Findroos et al. 2014) have the potential to 
help. There also remains the open question as to whether resourees to store the raw visibilities 
from SKA observations for later analysis will be available. 

In the near-term, we expeet to begin addressing these issues by extending previous Gravita¬ 
tional FEnsing Aeeuraey Testing (GREAT) ehallenges (e.g. Mandelbaum et al. 2014) to radio data 
with a radioGREAT ehallenge^. This will aim to better quantify the eurrent status of shape mea¬ 
surement teehniques for radio weak lensing, assess their behaviour aeross ehanging data parameters 
and spur new interest and developments in the field. 

2.4 Polarization and rotation velocities as indicators of intrinsic alignments 

One unique advantage of radio teleseopes for weak lensing is the polarization information 
whieh ean provide information on the intrinsie (unlensed) shapes of baekground galaxies. As 
deseribed in Brown & Battye (2011), the position angle of the integrated polarized emission from 
a baekground galaxy is unaffeeted by gravitational lensing. If the polarized emission (whieh is 

• http ://radiogreat .j b. man. ac .uk 
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polarized synchrotron emission sourced by the galaxy’s magnetic field) is also strongly correlated 
with the disk structure of the galaxy then measurements of the radio polarization position angle can 
be used as estimates of the galaxy’s intrinsic (unlensed) position angle. 

Such an approach could potentially have two key advantages over traditional weak lensing 
analyses. Firstly, the polarization technique can be used to effectively remove the primary astro- 
physical contaminant of weak lensing measurements - intrinsic galaxy alignments (Heavens et al. 
2000; Catelan et al. 2001; Hirata & Seljak 2004; Brown et al. 2002) - which are a severe worry for 
ongoing and future precision cosmology experiments based on weak lensing. Secondly, depending 
on the polarization properties of distant background disk galaxies, the polarization technique has 
the potential to reduce the effects of noise due to the intrinsic dispersion in galaxy shapes. 

The power of the polarization technique in practice will depend on two key observational 
parameters: the scatter in the relationship between the observed polarization position angle and the 
intrinsic structural position angle of the galaxy (Orms) and the number of galaxies for which one 
can obtain accurate polarization measurements (wpoi). These parameters depend on the details of 
the polarization properties of background galaxies (e.g. the mean polarization fraction, ITpoi) which 
are currently not well known. There are some existing measurements for a sample of local spiral 
galaxies (Stil et al. 2009) which suggest arms < 15° and ITpoi < 20% although the sample is small. 
Note that it may be possible to select sub-samples of the total galaxy population to have particular 
polarization properties. For example, one could imagine that selecting only galaxies with high 
fractional polarization would yield a galaxy sub-sample with highly ordered magnetic fields which 
would consequenfly have a very fighf correlation (low Orms) befween fhe polarization orienfafions 
and fhe infrinsic sfrucfural posifion angles of fhe galaxies. Of course, such a sub-sample would also 
have a very low surface number densify of galaxies. The polarization fechnique may Iherefore be 
heller suiled fo probing fhe shear power speclrum on large scales where high number densilies are 
nol required. 

A second novel idea lhal is well suiled lo radio observations is lo use rolalional velocity mea- 
suremenls lo provide information aboul Ihe inlrinsic shapes of galaxies. The idea, lirsl suggested by 
Blain (2002) and Morales (2006), is lo measure Ihe axis of rolalion of a disk galaxy and lo compare 
Ibis lo Ihe orienlalion of Ihe major axis of Ihe galaxy disk image. In Ihe absence of lensing, Ihese 
Iwo orienlalions should be perpendicular and measuring Ihe deparlure from perpendicularity di- 
reclly estimates Ihe shear field al Ihe galaxy’s position. Such an analysis would require commensal 
HI line observations which could in principle be done al no exlra cosl in terms of telescope time. 
The rotation velocity technique shares many of the characteristics of the polarization approach de¬ 
scribed above - in the limit of perfectly well-behaved disk galaxies, it is also free of shape noise 
and it can also be used to remove the contaminating effect of intrinsic galaxy alignments. In prac¬ 
tice, the degree to which the rotational velocity technique improves on standard methods will be 
dependent on observational parameters analogous to the ones for polarization discussed above. Re¬ 
cently, Huff et al. (2013) have proposed to extend this technique using the Tully-Fisher relation to 
calibrate the rotational velocity shear measurements and thus reduce the residual shape noise even 
further. 

Both the polarization technique and the rotation velocity approach are currently being tested 
as part of the SuperCLASS, CHILES and CHILES con Pol projects. They offer great promise for 
reducing the impact of shape noise and intrinsic alignments in radio weak lensing surveys. 
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3. Weak leasing with HI Intensity Mapping 


Weak gravitational leasing of high redshift 21cm emission has been the subject of several stud¬ 
ies focusing on the Epoch of Reionization (EoR) observations. In Zahn & Zaldarriaga (2006) and 
Metcalf & White (2009) it was shown that if the EoR is at redshift z ~ 8 or later, a large telescope 
like the SKA could measure the leasing power spectrum and constrain the standard cosmological 
parameters. The authors extended the Eourier-space quadratic estimator technique, which was first 
developed in Hu (2001) for CMB leasing observations to three dimensional observables, i.e. the 
21 cm intensity field I{d,z)- These sfudies did not consider 21 cm observations from redshifts after 
reionization when the average HI density in the Universe is much smaller. 

HI intensity mapping is a technique that has been proposed for measuring the distribution of 
HI gas before and during reionization (see Pritchard et al. in this volume) and measuring the BAO 
at redshifts of order unity (Chang et al. 2008, 2010; Seo et al. 2010; Masui et al. 2010; Ansari et al. 
2012; Battye et al. 2013; Chen 2012; Pober et al. 2013). In this technique, no attempt is made to 
detect individual objects. Instead the 21 cm emission is treated as a continuous three dimensional 
field. The angular resolution of the telescope need not be high enough to resolve individual galaxies 
which makes observations at high redshift possible with a reasonably sized telescope. 

A recent study (Pourtsidou & Metcalf 2014) extended the 21cm-intensity-mapping-lensing 
method further, taking into account the discreteness of galaxies, and investigated the possibility of 
measuring lensing at intermediate redshifts without resolving (in angular resolution) or even iden¬ 
tifying individual sources. Here we present an improved analysis of the signal-to-noise expected 
from an SKA-like telescope array, and demonstrate the performance of the SKA-Mid facility for 
SKAl and for SKA2. 

One of the first objectives of a 21cm lensing survey will be to measure the two-point statistics 
of the convergence field k(L,Zs) or, equivalently, the displacement field 5G(L,Zs), averaged over 
Zi. That is, 

Cf £dzP{k = L/r{z),z){W{z)f/E{z), (3.1) 

where W{z) = (r(Zs) — r(z))/r(zs), E(z) = H{z)/Hq. The expected error in the power spectrum 
(jS9Se averaging over L directions in a band of width AL is given by 


/^(jSQSe 


(2L+l)AL/sky 




(3.2) 


Here, Nl is the derived 5 G estimate reconstruction noise which involves the underlying dark 
matter power spectrum, the HI density Ohi(z) as well as the HI mass moments up to 4th order 
and, of course, the thermal noise of the telescope (see Pourtsidou & Metcalf 2014 for details). 
Eor the HI mass function we will first assume the most conservative scenario, i.e. a no evolution 
model using the results from the HIPASS survey (Zwaan et al. 2003), giving Hhi ~ 4.9 x 10^^ 
(Model A). We will also consider a more realistic evolution scenario from Zhang & Pen (2006) 
where = 4.9 x 10^^ (1 -|-z)^'^exp(—z/1.3). This is expressed through an evolution of the 

normalization of the Schechter function (Model B). 

Here we will concentrate on source redshift ~ 2. The SKAl-Mid array operates in the 
frequency range ~ 350 — 3050MHz divided in three bands (Dewdney et al. 2013). Our chosen 
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Zi corresponds to frequency = 473 MHz (Band 1). The power spectrum of the thermal noise 
of a SKA-like interferometer core array (using the uniform distribution approximation) is = 

—72—? where the system temperature T^ys = (40+ LlTskvlK with Tskv = B is 

the chosen frequency window, fobs the total observation time, Dtei the diameter (maximum baseline) 
of the array, £max(‘^) = 27rDtei/A is the highest multipole that can be measured by the array at 
frequency v (wavelength A), and /cover is the total collecting area of the telescopes Acoii divided 
by 7r(Dtei/2)^. For SKA-Mid we can consider a 2yr observation time, /sky ~ 0.7 (corresponding 
to ~ 30,900deg^ survey area), and we further choose B = 20MHz and AL = 36. Keeping these 
values constant, in Fig. 5 we present forecasts for the displacement field power spectrum for source 
redshifts = 2 and Zi = 3. For the Model A HI mass function, SKAl does not have the sensitivity to 
measure the signal. However, for Model B we see a great improvement in the signal-to-noise ratio, 
such that the lensing signal can be well measured by SKAl (Fig. 5, left panel). For SKA2, a precise 
characterization of the power spectrum is possible at multiple redshifts even in the conservation no¬ 
evolution scenario (Model A, Fig. 5, right panel). 




Figure 5: Left. Displacement field power spectrum for Zs = 2 and the corresponding measurement errors 
using the SKAl-Mid specifications and Model B for the HI mass function. Right. Displacement field power 
spectrum for Zj = 2 (solid black line) and z, = 3 (solid magenta line) and the corresponding measurement 
errors using the SKA2 specifications and Model A (no evolution) for the HI mass function. See the main 
text of Section 3 for details of the HI mass function adopted for Models A and B. 

This technique should enable us to measure the lensing power spectrum at source redshifts 
well beyond those accessible with more traditional weak lensing surveys based on the shearing of 
individual galaxy images. This is a promising approach to investigating the possible evolution of 
dark energy or modified gravify af high redshiff. In addition fo measuring fhe power specfrum, fhe 
SKA-Low array should be able fo map fhe dark matter disfribufion wifh high fidelify in a 5-by-5 
degree field using fhis technique (see Prifchard et al. in fhis volume for defails). 

4. Galaxy and cluster weak lensing 

In addifion fo probing large-scale sfrucfure, weak lensing can be used fo measure fhe disfribu- 
fion of mass (bofh dark and bright) around individual galaxies, as well as other structures such as 
galaxy pairs, groups and clusters, or voids (Mandelbaum et al. 2006; Gillis et al. 2013; Melchior 
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et al. 2014). Since the signal from an individual lens is too weak to be detectable, except for mas¬ 
sive clusters, such “galaxy-galaxy” measurements are usually stacked to produce an average mass 
distribution around a certain type of lens objects. On large, linear scales galaxy-galaxy lensing 
constrains galaxy bias (e.g. Hoekstra et al. 2002). Moreover, interpreted as a galaxy position-shear 
cross-correlation it can provide additional cosmological constraints and valuable calibration of as- 
trophysical systematics (Bernstein 2009; Joachimi & Bridle 2010) and Kirk et al, this volume). 
On smaller scales galaxy-galaxy lensing is sensitive to the radial matter density profile of galaxy 
haloes as well as the abundance and distribution of substructure in its outskirts, all as a function of 
galaxy properties such as type, stellar mass, or velocity dispersion (van Uitert et al. 2013). 

The SKA surveys will provide a unique angle to galaxy-galaxy weak lensing studies, with key 
contributions in the following areas: 

• Competitive or superior statistical errors: These errors are driven by the number density 
of both lens and source galaxies (particularly on small scales) and the sample variance due 
to a finite sky coverage (particularly on large scales). 

• Spectroscopic redshifts for new galaxy samples: To obtain unbiased measurements, a clean 
separation of lens and source samples is paramount. Spectroscopic redshifts for foreground 
objects allow for a clean and precise definition of lens samples and lens properties. Accurate 
redshifts for a subset of the source sample characterise the redshift distribution of background 
objects. 

• Novel lens galaxy sample definitions: The information from radio wavelengths can be used 
to define new lens samples e.g. in terms of their radio luminosity, AGN activity, or gas 
abundance. This allows for the measurement of the relation between these properties and 
matter halo characteristics, e.g. the HI gas mass to halo mass ratio as a function of redshift. 


To quantify this performance, we forecast expected errors on a one-parameter radial matter 
profile constrained via the average tangential shear measured in the range of 0.1 Mpc/h to 1 Mpc//i 
transverse separation from the lenses. We use a singular isothermal sphere (SIS) as a simple profile. 


Psis (f) 


iTiGr'^ 


(4.1) 


where r denotes three-dimensional distance, and where Oy is the velocity dispersion, which is the 
single free parameter. To facilitate comparison, we set Oy = 250 km/s for all calculations. The 
tangential shear signal for this profile is given by (assuming a spatially flat Universe, Bartelmann 
& Schneider 2001) 


(7+)W=2;r(^) ^ ( 1 - 


XL 


XL 


(4.2) 


where is the comoving distance to the lens, and the comoving distance to the source galaxy 
used to estimate the gravitational shear. 

In Fig. 6 we show the forecasted relative error on the SIS free parameter Oy for a selection of 
lens and source sample combinations. We have divided the lens samples into redshift bins of width 
0.1 and 0.2, respectively. It is possible to explore new terrain in galaxy-galaxy lensing already 
with early SKAl data by selecting lens samples using spectroscopic redshifts from an HI survey. 
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HI selected SKAlearly + DBS sources, 1000 deg^ '— ^ 
HI selected SKAlearly + DBS sources, 5000 deg^ '— ^ 
HI selected SKAl + Buclid sources, 5000 deg^ '— *- 
HI selected SKA2 + SKA2 sources, 5000 deg;: 

HI selected SKA2 + SKA2 sources, 30000 deg^ ' —*- 
Ha selected Buclid + SKA2 sources, 10000 deg^ '— *“ 
Ha selected Buclid + Buclid sources, 10000 deg'^ '—^ 



^lens 


Figure 6: Relative la error on the velocity dispersion, ff,,, of a singular isothermal sphere (SIS) prohle as 
a function of the redshift of the lensed objects, constrained via the average tangential shear signal in the 
transverse separation range G [0.1; 1] Mpc/h. The legend lists the lens sample selection, the source sample 
from which the shear signal is measured, and the assumed overlap in survey area. The error bars show the 
assumed width of the redshift bins. A 1 % constraint is indicated by the horizontal dashed line. 


Combining this with a source sample from an optical survey such as the Dark Energy Survey 
(DES), one can obtain percent-level constraints out to z ~ 0.6 for a few subsamples (e.g. selected 
on gas abundance). The full SKAl in combination with Euclid imaging for a source sample will 
decrease errors further by about a factor of five and extend the range where lenses can be studied 
to z = 0 . 8 . 

SKA2 will enable unique and powerful galaxy-galaxy lensing measurements without the need 
to rely on optical counterparts. As Eig. 6 shows, it will generate an unprecedented redshift baseline 
to study halo properties with a single facility. Out to z ~ 1.5 constraints are below the per mil 
level, which allow for tight constraints on complex models and/or the division into several sub¬ 
samples. SKA2 will also be on par with Euclid in measuring the galaxy-galaxy lensing signal of 
emission-line (mostly Ha) galaxies detected in the Euclid spectroscopic survey, providing inde¬ 
pendent corroboration for this high-redshift measurement. Note finally that the very deep SKA2 
surveys will allow the application of the novel techniques described in Section 2.4 as well as mag¬ 
nification analyses (e.g. Morrison et al. 2012) in a galaxy-galaxy lensing context. 


5. Concluding remarks 

The field of weak lensing with radio surveys is currently in its infancy but the prospect of 
SKA surveys coming online at the end of this decade presents a huge opportunity for opening a 
new wavelength window on the dark Universe. We have argued in this Chapter that the radio band 


15 











Weak tensing with the Square Kilometre Array 


M. L. Brown 


and the SKA in particular offer exciting prospects for enhancing and extending the reach of weak 
lensing studies beyond the limits of what is possible using their traditional optical approaches. 

The SKA surveys will probe the galaxy population to higher redshift where the lensing signal 
is larger and potentially easier to measure. There is huge potential to exploit the cross-correlation 
of lensing shear measurements obtained from overlapping optical and radio surveys which can be 
powerful in mitigating a number of systematic effects. 

Moreover, the radio band offers the exciting possibility of applying new techniques to measure 
weak lensing signals. One can make use radio polarization and/or rotational velocity information 
(from HI observations) to mitigate against intrinsic galaxy alignments which are the most worrying 
systematics for future precision weak lensing measurements. The opportunity to analyze HI Inten¬ 
sity Mapping surveys for the effects of weak lensing will open a new window on the dark Universe 
at redshifts that are well beyond the reach of traditional cosmic shear techniques. 

The SKA also offers a new route to performing galaxy-galaxy lensing where spectroscopic 
redshifts from HI observations and the extra diagnostics coming from the radio band will allow for 
the precise definition of new lens samples. 

As with other future large-scale survey instruments, weak lensing studies with the SKA will fa¬ 
cilitate the auto- and cross-correlation of galaxy clustering and cosmic shear measurements within 
the same experimental set up. Generally speaking, the joint exploitation of clustering and lensing is 
extremely useful as the two probes are highly complementary. Constraints from baryonic acoustic 
oscillations and/or redshift space distortions are “orthogonal” to weak lensing, especially in the 
presence of uncertainties in photometric redshifts and inaccurate knowledge of galaxy clustering 
bias (Zhan 2006; Abate et al. 2012; Laureijs et al. 2011; Amendola et al. 2013). 

In terms of constraining power for dark energy and/or modified gravity, galaxy clustering data 
measure the “dark fluid” equation of state at higher redshift than other standard candles like type la 
supernovae. Furthermore, clustering probes the evolution of matter fluctuations, which means that, 
through the Poisson equation, it is a proxy for the Newtonian potential <I>. In the case of modified 
gravity, clustering is sensitive to modifications occurring to the Newtonian gravitational constant. 
On the other hand, the deflection potential is proportional to the sum of the two metric potentials - 
which are equal in GR but not in more general gravity theories. Through this dependence, lensing 
effects represent a direct probe of non-standard gravity behaviours. As a consequence, the sensitiv¬ 
ity to beyond-GR growth parameters comes in its greatest part from weak lensing, which provides 
the only direct measurements of growth (without biasing), e.g. Weinberg et al. (2013). 

In order to realise the significant and exciting potential that radio weak lensing has to offer, 
a great deal of development work on analysis techniques, and demonstrating those techniques on 
real radio data from precursor surveys will be required. A key challenge will be to develop the 
field of galaxy shape esfimation for radio observations to the level of maturity that is currently 
enjoyed by the optical lensing community. Work has already started in this direction and will 
gain further momentum through the ongoing “radioGREAT” challenge during the coming year. 
At the same time, surveys with e-MERLIN, EOEAR and the JVEA are starting to produce the 
high-resolution, wide-held data suitable for demonstrating radio weak lensing techniques on real 
observations. The planned program of algorithm development work coupled with the successful 
application of the developed radio lensing tools to the SKA pathhnder and precursor surveys will 
position the community favourable to fully exploit the tremendous opportunity for radio weak 
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lensing that will be afforded by the eommissioning of the SKA at the end of this deeade. 
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